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ABSTRACT 

We present the results of continued monitoring of the quiescent neutron star low-mass X-ray binary 
XTE J1701— 462 with Chandra and Swift. A new Chandra observation from 2010 October extends 
our tracking of the neutron star surface temperature from ~800 days to ~1160 days since the end of 
an exceptionally luminous 19 month outburst. This observation indicates that the neutron star crust 
may still be slowly cooling toward thermal equilibrium with the core; another observation further 
into quiescence is needed to verify this. The shape of the overall cooling curve is consistent with 
that of a broken power law, although an exponential decay to a constant level cannot be excluded 
with the present data. To investigate possible low-level activity, we conducted a monitoring campaign 
of XTE J1701— 462 with Swift during 2010 April-October. Short-term flares — presumably arising 
from episodic low-level accretion — were observed up to a luminosity of '^1 x 10^^ erg s~^, ^^20 times 
higher than the normal quiescent level. We conclude that flares of this magnitude are not likely to 
have significantly affected the equilibrium temperature of the neutron star and are probably not able 
to have a measurable impact on the cooling curve. However, it is possible that brighter and longer 
periods of low-level activity have had an appreciable effect on the equilibrium temperature. 
Subject headings: accretion, accretion disks - stars: neutron - X-rays: binaries - X-rays: individual 
(XTE J1701-462) 



1. INTRODUCTION 

Observations of quiescent neutron star transients have 
the potential to constrain the internal properties of neu- 
tron stars. Some transients (so-called quasi-persistent 
transients) undergo years- or decades-long outbursts dur- 
ing which the neutron star crust is heated out of thermal 
equilibrium with the core due to nuclear reactions in- 
duced d eep in the crust (so-called deep crustal hea ting; 
see, e.g.. lBrown et aLlUQQafHaensel fc Zdunikll2008D : ob- 
serving the cooling of the surface after the source re- 
turns to quiescence can give information on the prop- 
erties of the crust (e.g., iRutledgc ct al. 2002). In ad- 
dition, the equilibrium crust and core temperatures in 
accreting neutron stars are set by the long-term time- 
averaged mass accretion rate and the efficiency of the 
cooling mechanisms (neutrino emission) at work in the 
neutron star interior; this cooling in turn depends on the 
properties of the matter inside the star (for an overview 
on n eutron star cooling, see, e.g., lYakovlev fc PethicS 
120041) . Measurement of the equilibrium surface tempera- 
ture can therefore yield information on the material deep 

^ Department of Physics, Massachusetts Institute of Technol- 
ogy, 77 Massachusetts Avenue, Cambridge, MA 02139, USA; 
joelkfOmit.edu 

^ MIT Kavh Institute for Astrophysics and Space Research, 
77 Massachusetts Avenue, Cambridge, MA 02139, USA 

^ Astronomical Institute "Anton Pannekoek," University of 
Amsterdam, Science Park 904, 1098 XH Amsterdam, The 
Netherlands 

* Institute of Astronomy, University of Cambridge, Madingley 
Road, Cambridge CB3 OHA, UK 

^ Department of Physics and Astronomy, Michigan State Uni- 
versity, 3250 Biomedical and Physical Sciences Building, East 
Lansing, MI 48824, USA 

^ Kapteyn Astronomical Institute, University of Groningen, 
P.O. Box 800, 9700 AV Groningen, The Netherlands 

INAF-Osservatorio Astronomico di Brera, Via E. Bianchi 46, 
1-23807 Merate (LC), Italy 



inside the neutron star. This, however, also requires good 
knowledge of the average mass accretion rate onto the 
star, and this rate is often poorly known. 

XTE J1701— 462 is a transient neutron star low-mass 
X-ray binary (NS-LMXB) that underwent an exception- 
ally luminous 19 month out burst in 2006-2007 ()Lin et al.l 
I2009bt iHoman et all I2O1O0 . We have monitored the 
source since it entered quiescence in 2007 August; in 
IPVid riksson et al. (2010, hereafter referred to as Paper I) 
we reported on Rossi X-ray Timing Explorer (RXTE) 
and Swift observations tracking the transition of the 
source from outburst to quiescence, and on Chandra and 
XMM-Newton observations covering the first ~800 days 
of the quiescent phase. The effective surface temperature 
of the neutron star — inferred from the thermal compo- 
nent of the quiescent spectra — exhibited a large decrease 
during the first ~200 days of quiescence, presumably due 
to the rapid initial cooling of the neutron star crust; such 
rapid cooling strongly indicates a highly conduc tive crust 
(|Shternin et al.|[20Ci7llBrown fc Cumminell2009l ). An ob- 
servation ~230 days into quiescence showed a large in- 
crease in both thermal and nonthermal flux from the 
source. A subsequent observation (~70 days later) sug- 
gested still somewhat elevated thermal flux; later obser- 
vations (^200 days after the initial increase) were con- 
sistent with slow cooling from the level preceding the 
increase. After ~800 days of quiescence, it was not clear 
whether the neutron star crust had already reached ther- 
mal equilibrium with the core or was still slowly cooling. 

XTE J1701-462 is one of only four NS-LMXBs for 
which the cooling of the neutron star crust has been 
monitored after long periods of accretion-induced heat- 
ing; the other three are KS 1731—260 (Wiina nds et al.l 
2001, 2002; Cackctt ct al. 2006, 2GlGa), MXB 16 59-2 9 
(IWiinands et al.u200a. .2004 ,Cackett et al..,200l [20081), 
and EXO 0748-676 (IDeeenaar et al.ll2009ll2QTTI) . In ad- 
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dition, a recent observation of the Terzan 5 globular clus- 
ter transient IGR J17480— 2446 indicates that the neu- 
tron star crust may have been heated significa ntly during 
a very bright ^2 month outburst (Degcnaar fc Wiinandsl 
I2011airbf ): further observations are needed to confirm this 
by observing cooling. As discussed in Paper I, the quies- 
cent behavior of XTE J1701— 462 has exhibited several 
noteworthy characteristics. Significantly higher effec- 
tive surface temperatures have been observed from XTE 
Jf70f— 462 than the other three sources; this may be 
partly due to a high core temperature in the neutron star. 
The cooling data from the first ~800 days of quiescence 
are well fitted — when excluding the observations affected 
by the temporary increase in temperature — with both an 
exponential decay to a constant level (e-folding time of 
~ 120^20 days) and with a broken p ower-law model. The 
latter is particularly relevant, since iBrown fc Cummingl 
(|2009| ) calculate that the general form of the cooling 
curve should be a broken power law leveling off to a 
constant at late times. However, the predicted break — 
mainly due to a change in heat capacity where the 
crust material transitions from a classical to a quantum 
crystal — is expected to occur a few hundred days into 
quiescence, whereas the observed break we reported in 
Paper I is at ~20-150 days, and may therefore have a 
different origin. In addition to the thermal spectral com- 
ponent, XTE J1701— 462 has shown a prominent and 
highly variable nonthermal component throughout the 
quiescent phase. Nonthermal components have been seen 
in the quiescent emission of many NS-LMXBs and are 
usually well fitted with a simple power law with a pho- 
ton index of 1~2. The origin of this nonthermal emission 
is poorly understood; pulsar shock emission or low-level 
accretion onto the neutron star surface o r magnetosphere 
have been suggested as explanations (jCampana et al.l 
119981 ). In Paper I, we speculated that the power- law 
component seen from XTE J1701— 462 and the large in- 
crease in luminosity more than 200 days into quiescence 
mentioned above are due to episodic low-level residual 
accretion. 

During 2010 April-October we conducted a monitor- 
ing program of XTE J1701— 462 using Swijt, with short 
(~3 ks) observations taking place once every two weeks. 
The goal was to study possible low-level activity during 
quiescence — as suggested by the temporary increase in 
luminosity seen during our Chandra and XMM-Newton 
monitoring — and investigate whether such activity can 
to some extent explain the high surface temperatures 
observed for XTE J1701-462. In 2010 October, we also 
obtained a new Chandra observation of the source to con- 
strain possible ongoing cooling of the neutron star crust. 
In this paper, we report on the results of our Swift mon- 
itoring campaign and recent Chandra observation. 

2. DATA REDUCTION AND RESULTS 

2.1. Cooling Analysis 

XTE J1701-462 was observ ed with Ch andra on 2010 
October 11-12. This ACIS-S ()Garmire et al. 2003) ob- 
servation (ObsID 11087) was performed in Timed Ex- 
posure mode, with the source located at the default S3 
aimpoint. The detector was operated in full-frame mode 
(with a frame time of 3.2 s), using the Very Faint teleme- 
try format. The observation had a live exposure time of 



56.60 ks; no periods of backgrou nd flaring were found. 
We analyzed the data with CIAO (jFruscione et al.|[20 06) , 
versioi i 4.2 (CALDB, ve r. 4.3.1), and with ACIS Ex- 
tract (|Broos et al.ll2010( ). version 2010-07-09. At the 
same time, we also re-analyzed all older Chandra obser- 
vations of the source in quiescence to ensure complete 
consistency. For information on the earlier observations, 
see Paper I. We mostly follow the same analysis proce- 
dure as described in Paper I and refer to the descrip- 
tion therein, mentioning here mainly any differences in 
the procedure. For each observation, background counts 
were extracted from an annulus with an inner radius of 
~5'.'3 (~10.8 pixels) and an outer radius of ~14'.'3 (~29.0 
pixels). All observations were checked to see whether 
applying a temperature-dependent charge transfer inef- 
ficiency correction to the data was appropriate. In all 
cases, the focal plane temperature stayed within half a 
degree of the nominal temperature of — 119.7°C through- 
out the observation, and therefore no correction was nec- 
essary. Two of our 11 Chandra observations consisted 
of two or more exposures spread over a few days. The 
spectra from the four subexposures of the seventh Chan- 
dra observation — taken over a period of ^3 days — were 
combined into a single spectrum using ACIS Extract, as 
were the response files for the exposures. This was done 
to facilitate fitting of the spectrum; the shortest ex- 
posure only had ~150 counts. As in Paper I, the Chan- 
dra spectra were grouped with the ACIS Extract tool 
ae_group_spectrum, with an average of ~25-30 counts 
per group. The tenth Chandra observation — consisting 
of two exposures performed over a ~27 hr period — was 
treated in the same way as the seventh one (note, how- 
ever, the discussion in Paper I of possible variability be- 
tween the individual exposures in both the seventh and 
tenth observations). 

We also re-analyzed the three XMM-Newton observa- 
tions previously reported on in Paper I with the latest 
software and calibration (SAS, ver. 10.0.2). The analy- 
sis proceeded in much the same way as described in our 
previous paper, and the description of identical analy- 
sis steps will not be repeated here. The XMM-Newton 
spectra were grouped with the specgroup task, requir- 
ing a minimum signal-to-noise ratio of 5, and limiting 
any oversampling of the intrinsic detector resolutions to 
less than a factor of 2.5. This resulted in an average 
number of counts per group of ^^35-40, except for the pn 
spectrum of the third XMM-Newton observation, which 
had an average of ~51 co unts p er group. 

Using XSPEC (|Arnaudlll996l) . version 12.6.0, we fitted 
all the Chandra and XMM-Newton spectra simultane- 
ously in the 0.5-10 keV band, except for the third XMM- 
Newton observation (XMM-3), which we fitted sepa- 
rately. As discussed in Paper I, this spectrum is very dif- 
ferent from the others, showing much higher nonthermal 
flux and a significant increase in thermal flux compared 
to the previous observation; XMM-3 will be discussed 
further below. We used the nsatm os neutron star atmo- 
sphere model (jHeinke et al.|[2006t see also a discussion 
of this model in Paper I) plus a simple power-law model 
(pegpwrlw), modified by photoelectric absorption. We 
used the TBnew abso rption modeP|, an updated version 
of the TBabs model (jWilms et al.lT2000D . with the vern 

* See |http : //pulsar ■ sternwarte .uni-erlangen. de/wllms/research/tbabs/ 1 
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cross sections ([Verner et al.lll996D and wilm abundances 
(jWilms et al.l2000r) . AlltheTBnew parameters except the 
equivalent hydrogen column density were kept fixed at 
their default values. The TBabs/TBnew mo del improves 
in va rious ways on older absorption models (jWilms et alj 
|2000( ) and is therefore to be preferred. In the nsatmos 
model, we fixed the neutron star mass at a va lue of 1.4 
Mq, the distance at 8.8 kpc ()Lin et al.ll2009al) . and the 
fraction of the neutron star surface emitting at 1. As 
in Paper I, we tied the absorption column and photon 
index between all observations (except for XMM-3). Al- 
though we have no reason to expect the photon index 
to necessarily have the same value for all the observa- 
tions, our spectra do not have enough counts to allow it 
to vary freely for each observation (see further discussion 
of this and its possible effects on our results in Paper I). 
We allowed the neutron star radius to float to its best-fit 
value initially (10.5 km), and then fixed it at this value 
when performing error scans for other free parameters 
(an error scan for the radius gives an uncertainty of ±2.5 
km). 

In Table [1] we show the (gravitationally redshifted) 
effective neutron star surface temperature for each in- 
dividual observation, as inferred from our spectral fits. 
Uncertainties in the neutron star radius, mass, and dis- 
tance give rise to a systematic uncertainty in the tem- 
peratures. This uncertainty is highly correlated between 
the different observations, and we do not include it in 
the error bars in Figures [T] and [2] and Tables [T] and [21 
Changes in the radius, mass, or distance systematically 
shift all the temperature values in more or less the same 
way (the higher the temperature, the larger the shift); in 
general, our conclusions about the cooling would not be 
significantly affected by such shifts. To assess the magni- 
tude of this uncertainty, we changed the values of these 
three parameters one at a time and re-fitted the spec- 
tra (allowing the absorption column and photon indices 
to float to new best-fit values). We consider a neutron 
star radius range of 8-15 km, a mass range of 1.1-2.0 
Mq, and a distance range of 7.5-10.1 kpc (bas ed on the 
1.3 kpc uncertainty in the distance estimate; iLin et aTl 
[2009a) . Changing the radius to 8 km (15 km) results in 
increases in the temperatures of ~6. 2-10.0 eV (decreases 
of ~9.0-13.5 eV). Changing the mass to 1.1 Mq (2.0 Mq) 
gives temperature increases of ~2. 5-4.1 eV (decreases of 
~5.5-8.6 eV). Finally, changing the distance to 7.5 kpc 
(10.1 kpc) yields decreases of ~6.3-9.2 eV (increases of 
~5.7-8.7 eV). Overall, we estimate that this gives rise 
to a combined systematic uncertainty of ^10 eV in the 
temperatures measured in recent observations. We stress 
that the parameter values derived from our cooling curve 
fits (see below) are not significantly affected by this (with 
the exception of the baseline temperature of the exponen- 
tial fit). This is further discussed and quantified in Paper 
I; there we also discuss uncertainty arising from the han- 
dling of the nonthermal component. We note here that 
all errors quoted in this paper correspond to la Gaussian 
(68.3%) confidence. 

In Table [1] we also show the 0.5-10 keV power-law 
flux and total luminosity, and the bolometric luminosity 
(all unabsorbed), as well as the fractional contribution of 
the power-law component to the total unabsorbed 0.5- 
10 keV fiux. We plot the luminosity, temperature, and 
power-law flux in Figure [1] The best-flt values for the 



equivalent hydrogen column density and the tied pho- 
ton index are (2.98 ±0.04) x 10^^ cm-^ and 1.96 ±0.20. 
The quality of the simultaneous fit is good: the reduced 
ixi) is 0.94 for 314 degrees of freedom (dof). The 
separate fit to XMM-3 is worse: xl = 1-30 (197 dof). 
In particular, this fit tends to underestimate the flux 
in the lowest energy bins, and the spectrum prefers a 
lower value for the absorption column. However, it is 
rather implausible that this observation would have sig- 
nificantly lower absorption than other observations dur- 
ing both quiescence and outburst (see below), and we 
therefore fix the column density at the best-fit value from 
the simultaneous fit (quoted above). The best-fit value 
of the XMM-3 photon index is 1.40 ± 0.08. We also tried 
fitting the XMM-3 spectrum with the power-law compo- 
nent replaced by the simpl or compTT Coniptonization 
models, but we were unable to improve on the fit. We 
note that it is conceivable that this spectrum is affected 
by the presence of elements heavier than hydrogen in the 
neutron star atmosphere. For an accretion rate above 
Mz 4 X 10-" Mq yr-i {8/Z){kTcs/0.1 keV)^/^, an 
element with atomic number Z is not expected to settle 
quickly enough out of the atmosph ere for it to main- 
tain a pure hydrog en composition (j Brown et al.l 119981 : 
iBildsten et"al1ll992[ ). The observed luminosity in XMM- 
3 corresponds to M ~ 2 x 10"^^ Mq yr~^ (assuming 
L — eM(?, with e = 0.2, and that ~80% of the emission 
arises from accretion) , implying that a presence of helium 
is possible if most or all of the additional flux — compared 
to the thermal emission of the previous observation — 
stems from accretion onto the neutron star surface. (We 
can also not exclude the possibility that some of the other 
observations are affected by heavier elements in the at- 
mosphere, although this is less likely; the fact that those 
spectra are in general well fitted by a pure hydrogen at- 
mosphere model indicates that this is probably not a 
large effect if present.) Another possible reason for the 
low quality of the XMM-3 fit is that thermal flux arising 
from ongoing accretion during this observation acts to 
distort the shape of the thermal component from that of 
a neutron star atmosphere spectrum. Finally, we note 
that it is co nceivable that cr o ss-cal ibration errors play 
a role here. iTsuiimoto et al.l (|2011f ) explore the cross- 
calibration between various X-ray detectors and derive 
lower absorption column values for XMM-Newton EPIC 
data, especially the pn detector, than for Chandra ACIS 
data. For XMM-3, the pn spectrum contributes the ma- 
jority of the residuals. However, even when allowing the 
column density to be free in the XMM-3 flt, the quality 
of the flt is still marginal {xt — 1-15 for 196 dof), and the 
change in the column de nsity (lower by c±14%) is some- 
what larger than seen bv lTsujimoto et al.l (pOlB (-11% 
for the pn and '^4%-6% for the MOS detectors). In ad- 
dition, the other two XMM-Newton observations agree 
well with the simultaneous fit and show no indication of 
the pn spectrum preferring a lower column density value 
than the MOS spectra. 

The temperature values reported here for the older 
Chandra and XMM-Newton observations are c±0.4-1.5 
eV lower than the corresponding values reported in Pa- 
per I. This is due to small differences in the data anal- 
ysis (value used for the neutron star radius; absorption 
model used; handling of XMM-3, CXO-7, and CXO-10; 
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Table 1 

Derived Parameter Values from Chandra and XMM-Newton Observations 



Observation'' 


t - to^ 
(days) 


(eV) (10- 


-13 Qj-g s 1 2 -J 


PL Fraction" 

(%) 


(10^3 erg s-i) 




CXO-1 


2.95 


163.1 ±3.5 


1.5 ± 1.1 


9±6 


16.6 ± 1.5 


15.7 ± 1.7 


CXO-2 


10.81 


158.2 ±2.5 


1.4 ±0.7 


9±4 


14.7 ± 1.0 


13.8 ± 1.2 


AiViiVi- _L 


16 24 




yj.o m u.o 




-i-O.u in u.o 




XMM-2 


49.49 


149.2 ± 1.3 


0.8 ± 0.3 


7±2 


11.6 ±0.4 


10.5 ± 0.5 


CXO-3 


174.33 


128.7 ±4.7 


4.9 ± 0.9 


47 ±6 


6.4 ±0.9 


9.6 ± 0.7 


XMM-3 


225.72 


157.8 ± 1.8 


15.2 ± 0.6 


53 ±2 


14.6 ±0.7 


26.5 ± 0.5 


CXO-4 


298.30 


135.1 ±2.0 


1.4 ±0.3 


17 ±3 


7.8 ±0.5 


7.6 ± 0.6 


CXO-5 


431.07 


125.5 ±3.1 


3.1 ± 0.5 


39 ±5 


5.8 ±0.6 


7.4 ±0.5 


CXO-6 


540.08 


125.0 ±1.5 


0.5 ± 0.1 


10 ±2 


5.7 ±0.3 


5.0 ±0.3 


CXO-7 


592.68 


128.4 ±2.3 


2.2 ± 0.4 


29 ±4 


6.4 ±0.5 


7.1 ± 0.4 


CXO-8 


652.62 


123.3 ±2.2 


1.8 ±0.3 


28 ±4 


5.4 ±0.4 


5.9 ± 0.3 


CXO-9 


705.38 


123.0 ±2.0 


1.5 ± 0.3 


25 ±3 


5.4 ±0.3 


5.6 ±0.3 


CXO-10 


795.63 


123.7 ±1.7 


1.1 ± 0.2 


19 ±3 


5.5 ±0.3 


5.3 ±0.3 


CXO-11 


1158.84 


120.6 ±1.7 


0.8 ±0.2 


16 ±3 


5.0 ±0.3 


4.6 ±0.3 


Note. — All numbers were derived assuming 


; a neutron star radius of 10.53 km, 


a mass of 1.4 Mq 


, and a distance 



of 8.8 kpc. Errors quoted are at the la Gaussian (68.3%) confidence level. 
^ See IPridriksson et al.l 120101) for technical information on all observations except the latest one. 
^ Time of mid observation; to is MJD 54321.95. 

Effective surface temperature of the neutron star as seen by an observer at infinity. 

Unabsorbed 0.5-10 keV flux of the power-lawr component. 
° Fractional contribution of the power-law component to the total unabsorbed 0.5-10 keV flux. 
' Unabsorbed bolometric luminosity of the thermal (nsatmos) component, 
s Unabsorbed total luminosity in the 0.5-10 keV band. 



Table 2 

Best-fit Cooling Curve Parameters 





Exponential Decay Fit 






Broken Power-law Fit 




Excluded^ 


(days) 


(eV) 


kT"= 
(oV) 


Xt (dot)" 




72' 


(days) 


A" 

(eV) 


xt (dof)° 




230 ± 46 


123.4 ± 0.9 
121.9 ± 1.5 


36.9 ± 1.7 
35.8 ± 1.4 


1.07 (9) 
1.66 (10) 


0.030 ± 0.013 
0.030 ± 0.012 


0.069 ± 0.004 
0.069 ± 0.004 


00 + 24 

■^8-12 

™-13 


168.8 ± 5.7 
168.8 ± 5.5 


0.88 (8) 
1.13 (9) 


XMM-3 and CXO-4 
XMM-3 



Note. — Errors quoted are at the Icr Gaussian (68.3%) confidence level. 
^ Best-fit e-folding time of the decay. 
^ Best-fit constant offset to the decay. 

Best-fit decay amplitude. 
^ Reduced for the fit and number of degrees of freedom. 
° Best-fit pre-break power-law slope. 
^ Best-fit post-break power-law slope. 
^ Best-fit break time between power laws. 

Best-fit pre-break normalization coefficient (i.e., temperature at i — io — 1 day). 
^ Observations excluded from the fit. 



more recent calibration) and the effect of the new obser- 
vation on the simuhaneous fit. The best- fit column den- 
sity value quoted above is significantly higher than the 
one reported in Paper I, primarily due to the different 
abundance table used here. We note that when fitting 
with the phabs absorption model with angr abundances 
and heme cross sections, as we did in Paper I, we get 
an absorption colum n of 2.0 x lOf^cm"^, the same as 
the value derived bv lLin et all ([2009b) from RXTE and 
Swift observations during the 2006-2007 outburst (and 
very close to the value derived in Paper I). 

The new Chandra observation — made ~1160 days into 
quiescence — has an effective neutron star surface tem- 
perature (as seen by an observer at infinity) kT^ = 
120.6 ± 1.7 eV, compared to 123.7 ± 1.7 eV in the previ- 
ous Chandra observation a year earlier. These two tem- 
perature values differ at the ~1.3cr (i.e., ^^80%) confi- 
dence level; using the weighted temperature average for 
the three Chandra observations preceding the latest one, 
123.4 ± 1.1 eV, gives similar significance for a change 



in the temperature. The data are therefore more con- 
sistent with continued cooling than with the tempera- 
ture having reached a fixed equilibrium value, but given 
the size of the error bars it is not possible to draw a 
firm conclusion. We also note that uncertainties asso- 
ciated with the handling of the nonthermal component 
can possibly give rise to small shifts in the tempera- 
tures (see discussion in Paper I). The current temper- 
ature corresponds to a bolometric thermal luminosity of 
2±5.0 X 10'^^ erg s~^. The power-law flux of the new ob- 
servation is among the lowest seen during quiescence and 
represents ~16% of the total 0.5-10 keV unabsorbed fiux. 
The 0.5-10 keV unabsorbed luminosity of the new obser- 
vation (~4.6 X 10^'^ erg s~^) is the lowest observed so far 
from the source. 

To estimate the end time of the outburst, ipj we fol- 
lowed the same procedure as in Paper I. We fitted the lu- 
minosity measured in the three Swift observations made 
during the final decay of the outburst with a simple ex- 
ponential decay function, and did the same for the first 
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Figure 1. Spectral evolution of XTE J1701— 462 since the end 
of the 2006-2007 outburst. Top panel: unabsorbed luminosity in 
the 0.5-10 keV band. We do not show a Scr Swift upper limit to 
the luminosity (observation 16; shown in Figure [Sj, which is vir- 
tually coincident with the last Chandra data point. Middle panel: 
effective neutron star surface temperature. The solid line is the 
best-fit exponential decay to a constant level, excluding the sixth 
and seventh data points. The dashed line is a similar fit including 
the seventh data point. Bottom panel: unabsorbed 0.5—10 keV 
power-law flux. 

three Chandra and XMM-Newton observations in quies- 
cence; we then define to as the intersection of the two 
curves. We did not re-extract the Swift spectra from Pa- 
per I but did re-fit them. The first two spectra were 
grouped with a minimum of 25 counts per bin and fitted 
with a simple absorbed power law (TBnew*pegpwrlw), 
fixing the absorption column at the value derived from 
our observations in quiescence (quoted above) . The third 
spectrum only had 52 counts, and we therefore fitted 
the unbinned spectrum with the M^-statistic in XSPEC. 
In this case, we added an nsatmos component to the 
spectral model, since the luminosity of this observation 
is low enough for thermal radiation from the neutron 
star surface to make a significant contribution to the 
flux (we derive a 0.5-10 keV unabsorbed luminosity of 
(3.8±0.7) X 10^* erg s'^). We fixed the effective neutron 
star surface temperature at a value of log(rcff/ii') — 6.4 
(i.e., a gravitationally redshifted value of kT^ w 166 eV, 
slightly higher than the value measured in the first quies- 
cent observation) and the radius at 10.5 km (as derived 
from our quiescent spectra). Due to the low number of 
counts in the spectrum of the third observation, we also 
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Figure 2. Effective neutron star surface temperature during qui- 
escence, with best-fit cooling curves shown. The solid curve is a 
broken power law fitted to all data points except the sixth and 
seventh ones; including data point 7 in the fit yields a virtually 
identical curve. The dashed curve is the best-fit exponential de- 
cay curve with a constant offset, where points 6 and 7 have been 
excluded from the fit (shown as a solid curve in Figure [TJ. 

in that case fixed the photon index at a value of 2 (similar 
to the value measured in the previous Swift observation 
and the combined value in our quiescent fit). The result- 
ing value for to is MJD 54321.95, i.e., ~1.4 days after 
the final Swift observation and ~3.0 days before the first 
Chandra observation. This value for to is 0.18 days earlier 
than the value derived in Paper I. The difference arises 
mostly from a somewhat lower luminosity value derived 
here for the last of the three Swift observations; in Paper 
I we did not include an nsatmos component in the model 
for that spectrum. 

As in Paper I, we fitted our derived temperatures with 
two models: an exponential decay function with a con- 
stant offset, T^it) = r'exp[-(i - to)/r] + T^q, and a 
broken power-law model (as a function of t — to). In both 
cases, we kept to fixed at the val ue quoted above. Th e 
fits were performed with Sherpa (jFreeman et al.|[2001f) : 
errors were estimated using the confidence method0 In 
all our fits we exclude XMM-3; as mentioned above, this 
observation has a much higher luminosity (arising from 
increases in both thermal and nonthermal flux) than the 
preceding and subsequent observations, likely due to an 
accretion event. We perform fits both with and without 
the subsequent observation, CXO-4. Although this ob- 
servation shows some evidence of increased surface tem- 
perature, it has neither anomalously high luminosity nor 
nonthermal flux. It is not clear whether this possibly 
elevated temperature is somehow due to increased accre- 
tion (before or during the observation), nor is it clear 
whether it is related to the elevated luminosity observed 
in XMM-3. It may possibly be a statistical fluctuation, 
or it may be due to our limited ability to separate the 
contributions of the two spectral components (see fur- 
ther discussion in Paper I as well as the discussion of 
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Table 3 

Swift XRT Observations of XTE J1701-462 in 
Quiescence 



No. 


Start Date'' 


ObsID 


Exp. Time"^ (ks) 


1 


Apr 2 


00090523001 


2.90 


2 


Apr 16 


00090523002 


3.43 


3 


Apr 30 


00090523003 


3.36 


4 


May 14 


00090523004 


3.06 


5 


May 27 


00090523005 


2.81 


6 


Jun 11 


00090523006 


2.85 


7 


Jun 25 


00090523007 


2.79 


8 


Jul 9 


00090523008 


3.05 


9 


Jul 23 


00090523009 


3.20 


10 


Jul 27 


00031776001 


3.00 


11 


Jul 31 


00031776002 


3.16 


12 


Aug 6 


00090523010 


2.66 


13 


Aug 20 


00090523011 


3.26 


14 


Sep 3 


00090523012 


1.56 


15 


Sep 17 


00090523013 


3.09 


16 


Oct 1 


00090523014 


3.08 


17 


Oct 27 


00090523015 


2.11 



All dates are in 2010. 
^ Good live exposure time. 

possible effects from low- level accretion in Section 1X4)) . 
The results of our cooling fits are shown in Table [2] and 
best-fit cooling curves are plotted in Figures [T] and [2j We 
note that the value used for to cannot affect the e-folding 
time or the equilibrium temperature of the exponential 
fits. It does affect the decay amplitude, but within the 
small uncertainty in to that effect is negligible (of the 
order of half an eV). However, the break time, and espe- 
cially the pre-break slope, of the broken power-law fits 
are quite sensitive to the value of to', the post-break slope 
is practically unchanged within the allowed range for to. 
We performed broken power-law fits (both XMM-3 and 
CXO-4 excluded) with shifted within the allowed range 
set by the last/first observation in outburst/quiescence. 
Shifting to to two days later (one day earlier) gives best- 
fit break times of ~30 (~42) days and pre-break slopes of 
~0.019 (~0.034). The effects of changes in to on broken 
power-law fits are further discussed in Paper I. 

Both the exponential and the broken power-law model 
give acceptable fits to the temperatures when both 
XMM-3 and CXO-4 are excluded: xl = 1-07 (9 dof) 
for the exponential and xl — 0-88 (8 dof) for the broken 
power law. However, an ongoing decrease in the tem- 
perature is inconsistent with the exponential fit. When 
CXO-4 is included, the broken power-law fit is superior: 
xl = 1.13 (9 dof), compared to xl = 1-66 (10 dof) for 
the exponential fit. The exponential fit with CXO-4 in- 
cluded (dashed line in Figure [1]) matches the later data 
points (t — to > 250 days) quite well but fits the first five 
data points badly. In contrast to the broken power-law 
fits, a single unbroken power law gives poor fits, with 
xl = 2.07 (10 dof) when both XMM-3 and CXO-4 are 
excluded; the best-fit slope is 0.059±0.002. We note that 
we cannot exclude the possibility that our inferred tem- 
peratures (and thereby the cooling curve fits) are affe cted 
by residual accretion; we discuss this in Section 15^ and 
argue that the effects of this, if present, are likely small. 

2.2. Swift Monitoring Program 

XTE J1701-462 was observed 17 times w ith the Swift 
X-ray Telescope (XRT: "B ^ws et al.|[2005l) during 2010 
April-October (see Table [3] for details on each observa- 
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Figure 3. Swift XRT light curve of XTE J1701-462 during 2010 
April-October (filled squares); the arrow represents a 3cr upper 
limit. The open circle indicates a contemporaneous Chandra obser- 
vation whose corresponding Swift XRT count rate was calculated 
from its spectrum; this count rate level is also indicated by the 
dashed line. 



tion). Fifteen of these were regularly scheduled obser- 
vations made approximately every two weeks as part of 
our monitoring program of the source; two additional 
Director's Discretionary Time (DDT) observations (nos. 
10 and 11 in Table [3]) were made in late July. All the 
observations were performed with the detector in Pho- 
ton Counting mode. We analyzed the data with HEA- 
SOFT, version 6.10. Using the xrtpipeline task, we 
processed and screened (with the default settings) the 
Level 1 raw event files and created vignetting-corrected 
exposure maps for each observation. We then extracted 
counts (and spectra) in the 0.5-10 keV range from the 
Level 2 screened event files using Xselect; source counts 
were extracted from a circle of radius 30" and back- 
ground counts were in each case extracted from a nearby 
circle of radius 250". For each observation, we used XIM- 
AGE to integrate the exposure map over both the source 
and background extraction circles, and used the ratio of 
the integrated exposures as our background scaling fac- 
tor, thereby taking into account the effects of vignetting 
and CCD bad pixels and hot columns. We also used the 
xrtmkarf task (with an exposure map) to calculate a 
correction factor for each observation to correct the net 
(i.e., background-subtracted) count rate for losses due to 
the finite size of the extraction region compared to the 
point-spread function (PSF), as well as for losses due to 
vignetting and CCD bad pixels and hot columns. The 
count rates we quote therefore correspond to an on-axis 
extraction region of infinite size. 

In Figure |Sl we show the Swift light curve of XTE 
J1701— 462. For one observation, in which the source was 
not detected, we show a 3(t upper limit. We also show 
a data point representing the latest Chandra observation 
(CXO-11); we used XSPEC to calculate the Swift XRT 
count rate corresponding to the best-fit Chandra spec- 
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trum. This CXO-11 count rate level is also indicated 
with a gray dashed line; we note that the spectrum from 
the previous Chandra observation — which was performed 
around 800 days into quiescence and had slightly higher 
thermal and nonthermal flux — corresponds to a ^40% 
higher count rate. Most of the Swift observations show 
count rates consistent with the quiescent level observed 
in CXO-11. However, the first Swift observation shows a 
clearly elevated count rate (a factor of ~6 increase com- 
pared to the expected quiescent level) and the ninth ob- 
servation shows a much stronger (factor of ^35) increase. 
The two following DDT observations, performed ~4 and 
~8 days later, show a rapidly decaying flux. An observa- 
tion performed two weeks after the brightest (i.e., ninth) 
observation shows a weak indication of an elevated count 
rate, but is consistent with the baseline quiescent level. 

For the three brightest Swift observations (nos. 1, 
9, and 10), we also performed spectral fitting (using 
XSPEC) in the 0.5-10 keV band. We created ancil- 
lary response files using the xrtmkarf task and used a 
standard redistribution matrix file from the calibration 
database. We used the corrected background scaling fac- 
tor described above. For observations 1 and 10, which 
only had 28 and 50 counts, respectively, we fitted the 
unbinned spectra with the T^-statistic in XSPEC. We 
binned the spectrum for observation 9 into groups with 
a minimum of 26 counts. We used the same spectral 
model as for the Chandra and XMM-Newton spectra, 
fixing all the parameters except the power-law normal- 
ization (but including the effective surface temperature) 
at the best-fit values for CXO-11 as derived from the 
simultaneous Chandra and XMM-Newton fit described 
above. For the brightest Swift observation (no. 9), we 
also allowed the photon index to vary; this gave a value 
of 2.11 ± 0.19, similar to the combined value from our 
Chandra and XMM-Newton fit. For all other Swift obser- 
vations (i.e., except the three brightest ones), we estimate 
a luminosity from the observed count rate, assuming a 
fixed thermal component equal to that seen in CXO-11 
and a varying power-law component with a photon in- 
dex value equal to that obtained from our simultaneous 
fit (i.e., we subtract the count rate corresponding to the 
absorbed CXO-11 thermal component — calculated with 
XSPEC — from the observed count rate and attribute the 
rest to an absorbed power-law component with photon 
index 1.96; we then convert this power-law count rate to 
an unabsorbed flux and add to the CXO-11 unabsorbed 
thermal flux). We show the Swift luminosities in Fig- 
ure [2 The peak luminosity observed in the 2010 July 
flare is ~1.0 x 10^^ erg s~^, i.e., a factor of ^20 higher 
than the expected quiescent level (as seen in our two 
most recent Chandra observations). The luminosity of 
Swift observations 9-13 (i.e., the initial peak detection 
of the flare and the four subsequent observations) can be 
fltted {xt = 0.13 for 3 dof) with an exponential decay 
down to a constant level of ~4.6 x lO"^'^ erg s~^, equal 
to that in the latest Chandra observation (which took 
place roughly 2 months later); the best-fit e- folding time 
is ~2.8 days. 

The limited number of counts in our Swift spectra make 
it difficult to place constraints on the spectral composi- 
tion and determine to what extent the increases in lu- 
minosity are due to the thermal/nonthcrmal component. 



Performing an error scan with the temperature free (in 
addition to the power-law normalization and slope) for 
the spectrum of the brightest Swift observation (no. 9, 
with ~180 counts) yields a 90% confidence upper limit 
on the effective temperature of kT^ c± 199 eV (and no 
lower limit). Fixing the temperature at that value and 
re-fltting the spectrum results in a fractional contribu- 
tion of the thermal flux to the total unabsorbed 0.5-10 
keV flux of ~27%. Instead flxing the power-law normal- 
ization at its 90% confldence lower limit and re-fltting 
gives a very similar result: a thermal fraction of ~26%. 
The data therefore seem to imply a 90% upper limit on 
the thermal contribution to the flux of ~30%. With the 
thermal flux around this upper limit, the flt prefers a 
lower photon index for the power law, ~1.3-1.4, similar 
to the value for XMM-3. 

Since none of the other Swift spectra have more than 
50 counts, we opt to explore this issue further with a 
simple hardness ratio analysis, rather than try to place 
fltting constraints on spectra with very few counts. We 
deflne our hardness ratio as the net counts in the 3-10 
keV band divided by those in the 0.5-3 keV band. This 
ratio is sensitive to the relative contributions of the two 
components, since the contribution of the thermal com- 
ponent above 3 keV is small, but a power-law component 
contributes signiflcantly there. In Figure [H we plot this 
hardness ratio against the (PSF-corrected) net count rate 
for the four brightest Swift observations (flUed squares; 
number of observation indicated next to the symbol). 
We also use XSPEC to calculate curves based on the 
TBnew* (nsatmos+pegpwrlaw) model and plot them in 
Figure HJ these curves depict various combinations of 
thermal and nonthermal flux in this hardness-intensity 
diagram (HID). In all cases, we assume the same values 
for the absorption column, distance, and neutron star 
mass and radius as in our simultaneous Chandra and 
XMM-Newton flt. The bottommost solid curve corre- 
sponds to pure thermal spectra; the other three solid 
curves represent spectra with flxed ratios of thermal and 
nonthermal flux (i.e., along these curves both the ther- 
mal and nonthermal components are varying, but in 
such a way that the relative contributions of the com- 
ponents to the total unabsorbed 0.5-10 keV flux stays 
constant) ranging from 25% power-law contribution (bot- 
tom) to 75% contribution (top). We calculate these four 
curves for temperatures ranging from that of CXO-11 
[kT^ = 120.6 eV) to the maximum temperature for 
which the nsatmos model is valid (log(roff/-?^) = 6.5; 
kT^ K. 210 eV); the 75% curve extends beyond the plot. 
For the three curves with a power-law contribution, we 
assume a photon index of F = 1.7; this is a compro- 
mise between the index seen for XMM-3 (1.40) and the 
tied index for the other observations in our simultaneous 
Chandra and XMM-Newton flt (1.96; we note that the 
best-flt index value for the penultimate Swift observa- 
tion during outburst at ~6 x 10'^^ erg s~^ was also very 
similar, 1.94 ± 0.10). We also plot curves correspond- 
ing to a flxed thermal flux plus a varying power-law flux 
for two different temperatures: 120.6 eV (as measured for 
CXO-11; long-dashed curves) and 157.8 eV (as measured 
for XMM-3; short-dashed curves). In the former case we 
do this for two photon index values: 1.7 (black curve) 
and 2 (gray curve); in the latter case for three values: 
1.7 (black), 1.4 (gray), and 2 (also gray). The intersec- 
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Figure 4. Hardness— intensity diagram for the four brightest Swift observations (number of observation indicated next to data point) . Also 
shown are the expected locations of the CXO-11 and XMM-3 spectra, as well as curves corresponding to spectra with various combinations 
of thermal (nsatmos) and nonthermal (power-law) contributions. Along the solid curves, the fractional contribution (indicated next to 
curve) of the power-law component to the 0.5-10 kcV unabsorbcd flux is constant while both spectral components vary (a photon index 
r = 1.7 is assumed). Along the dashed curves, the thermal flux is constant (with two different temperature values considered), but the 
power-law flux varies (the assumed value of the photon index is indicated next to each curve). Further details arc given in the text. 



tions of the solid curves and the F = 1.7 dashed curves 
indicate where along those dashed curves the particular 
power-law fraction has been reached. Finally, for com- 
parison wo also show the locations of the CXO-11 (open 
circle) and XMM-3 (cross) spectra in this plot. As can 
be seen from the dashed curves, the value assumed for 
the photon index has a significant impact on the curves. 
Changing the photon index would also shift the solid 
curves (except the purely thermal one) up (lower index) 
or down (higher index); the higher the power-law frac- 
tion, the larger the shift. The location of the XMM-3 
data point gives an indication of the magnittide of the 
shift for the 50% curve, since the XMM-3 spectrum has 
a power-law fraction of ~53% and a photon index of 1.40. 

The HID indicates that 5'wi/t observations 1, 10, and 11 
(especially no. 1) are not consistent with the count rate 
increase being purely due to an increase in the power-law 
component on top of the underlying baseline quiescent 
(CXO-11) thermal flux; there has to have been an in- 
crease in the thermal flux as well, given the location of 
these three data points to the right of the long-dashed 
curves. This is in agreement with the XMM-3 spectrum. 
However, the location of observation 9 in the diagram 
indicates that the power-law component dominates in 
that case, and this is in qualitative agreement with the 
spectral constraints on the thermal flux in this observa- 
tion (as discussed above). This implies that the increase 
in the thermal emission during a flare cannot keep pace 
with that of the nonthermal emission up to fluxes that 
high. Observations 9, 10, and 11 are consistent with a 
flare decay which starts with a decrease in power-law flux 



(F ^ 1.5 2) on top of a fixed thermal component with 
a temperature similar to that of XMM-3 (thus initially 
following a trajectory downward with a shape similar to 
the short-dashed curves); after the power-law fraction 
reaches ~5%-25%, the thermal flux decreases alongside 
the nonthermal flux (the trajectory then curving to the 
left and becoming closer to horizontal) until the base- 
line (CXO-11) temperature has been reached (i.e., the 
trajectory intersecting with the long-dashed curve). Ob- 
servation 1, although not part of the same flare, is also 
consistent with this. The location of XMM-3 perhaps 
suggests that the photon index of the power-law compo- 
nent in that case (1.40) was lower than in the Swift flare; 
however, it is important to keep in mind that the loca- 
tions of the two Chandra and XMM-Newton data points 
in the HID (as well as the plotted curves, of course) are 
model dependent, in contrast to the Swift data points, 
and this can affect comparisons. We also emphasize that 
any discussion of possible trajectories of flares in the HID 
based on the currently available data is obviously spec- 
ulative. 

3. DISCUSSION 

3.1. Cooling Curve 

We have presented a Chandra observation of the NS- 
LMXB XTE J1701-462, made in 2010 October, which 

extends our coverage of the source's current quiescent 
phase from ~800 days to ~1160 days after the end 
of the 2006-2007 outburst. This new observation sug- 
gests that the effective neutron star surface tempera- 
ture has decreased compared to the preceding Chan- 
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dra observation(s), implying that the neutron star crust 
may not have reached thermal equilibrium with the core 
yet. An additional observation — e.g., ~2000 days into 
quiescence — is needed to conclusively determine whether 
cooling is still ongoing, and if so, to constrain its rate. 
The new Chandra observation is consistent with the 
broken power-law fit to the cooling curve presented in 
Paper I. We cannot exclude the possibility that the 
cooling has followed an exponential decay down to a 
constant temperature level, although ongoing cooling 
would make this somewha t unlikely. The simulations of 
IBrown fc CummineJ ()2009( ) indicate that the form of the 
cooling curve should approximately be that of a broken 
power law leveling off to a constant at late times; how- 
ever, this break is predicted a few hundred days into qui- 
escence and is difficult to reconcile with the much earlier 
break in our power-law fits (see further discussion of this 
in Paper I). The fits presented in this paper indicate a 
break in the time range '^25-80 days post-outburst; if 
the temperature leveled off from the power-law behavior 
after ~200-300 days (which is somewhat unlikely given 
the new Chandra observation) , then the break could con- 
ceivably have been a few tens of days later. The range 
of break times derived here is similar to and consistent 
with that derived in Pap er I. The break predicted by 
IBrown fc Cummind (|2009[ ) is mainly due to a change in 
the heat capacity of the neutron star crust where the ma- 
terial transitions from a classical to a quantum crystal. 
A possible alternative explanation for the observed break 
in the XTE J1701— 462 cooling curve is the existence of a 
strong nuclear heating source in the crust not considered 
in the cooling curve sim ulations; in Paper I w e speculated 
that the fusion of ^^O (jHorowitz et al.ll2008D might be a 
possible candidate. As discussed in Paper I, the rapid 
initial cooling during the first ^200 days of quiescence 
indicates a neutron star crust with high thermal con- 
ductivity, and possibly suggests low-impurity material 
(IShternin et al.lf2007t lB rown fc Cummind[2009l) . Fitting 
theoretical models to the cooling curve is necessary to 
explore in quantitative detail the implications of the ob- 
served cooling on the internal properties of the neutron 
star; this is beyond the scope of this paper. We note 
again that it is possible that our results are affected to 
some extent by low-level accretion contributing to the 
observed thermal emission we use to infer the tempera- 
tures; we discuss this in Section [33] and conclude that 
such effects are unlikely to be large. 

It is instructive to compare the observed cooling 
curve of XTE J1701-462 to those of the other cool- 
ing quasi-persistent transients. The cooling curve of 
MXB 1659-29 can be fit ted weU with the model of 
IBrown fc Gumming! (j2009t ): this fit implies a break at 
~300-400 days post-outburst. However, a simple expo- 
nential decay to a constant level als o provides a good 
fit to the data (jCackett et al.l I2008D . Both fits indi- 
cate that the temperature reached a roughly constant 
level '^lOOO-lSOO days into quiescence. We note that 
the implied power-law slopes for MXB 1659—29 are 
much steeper than those of XTE J1701-462. We fitted 
the temperatures of the first four observations of MXB 
1659—29 (those made before the cooling leveled off) with 
a broken power l a w, us ing the temperature values given 
in iCackett et al.l ()2008l ) ; this gave pre- and post-break 
slopes of ~0.15 and ~0.55 (compared to '^0.03 and ~0.07 



for XTE J1701— 462). However, the data are very sparse 
for the early part of quiescence — the first observation was 
made ~31-38 days after the end of the outburst and the 
second observation did not take place until ~400 days 
into quiescence — and assumptions about the behavior of 
the source during that period should therefore be re- 
garded with caution. In contrast to MXB 1659—29, the 
coo hng curve of KS 1731— 260 is not well fitted by either 
the IBrown fc Cumminsj (|2009[ ) model or an exponential 
decay, whereas a single unbroken power law provides a 
good fit, and the source seems to still be slowly cooling 
more than 3000 days into quiesc ence: the slope of th e 
best-fit power law is 0.125 ±0.007 (jCackett et al.ll2010aD . 
However, the first observation of KS 1731—260 in quies- 
cence did not take place until ~48-65 days after the end 
of the outburst, and therefore an early break in the power 
law similar to the one indicated by the XTE J1701— 462 
data cannot be excluded. The four Chandra observa- 
tions of EXO 0748-676 made during the first ~600 days 
of quiescence can be fitted both with an exponential de- 
cay to a constant level and with a single unbroken power 
law (jPegenaar et al.|[2011l : the authors also analyze Swift 
observations of the source, which are less constraining 
but give results consistent with those from the Chandra 
observations). The best-fit power-law slope is 0.03±0.01. 
An early break cannot be excluded for EXO 0748—676 
either, since the first Chandra observation is only con- 
strained to have taken place sometime during the first 
~ 15-60 days of quiescence. 

In connection with the discussion above about the dif- 
ferent power-law slopes obs erved for the four s o urces , 
it is worthwhile to note that IBrown fc Cummind ()2009f ) 
show that the initial slope of the broken power-law curve 
gives a direct measure of the inward flux near the top 
of the crust during outburst. They also point out that 
the inferred early-time slopes for KS 1731—260 and MXB 
1659—29 imply a flux well in excess of that available from 
electron captures in the outermost layers of the crust. In 
contrast, the pre-break slope for XTE J1701— 462 (cou- 
pled with the mass accretion rate inferred from the ob- 
served outburst luminosity) implies a flux which is con- 
sistent with the energy available (see Paper I) . A possible 
explanation for the apparent inconsistency in the cases 
of KS 1731-260 and MXB 1659-29 might be that an 
early break such as that observed in XTE J 170 1—462 
took place in those sources as well. 

3.2. Low-level Activity 

We have also presented results from a Swift XRT 
monitoring campaign of XTE J1701— 462 consisting of 
17 short observations made during 2010 April-October. 
This campaign has clearly established that the increase in 
flux seen in our third XMM-Newton observation was not 
an isolated event, and that signiflcant temporary eleva- 
tions of the source flux are repeatedly taking place during 
quiescence (in the context of XTE J1701— 462, we refer 
to any period outside of extended bright outbursts as 
quiescence). In 2010 July, we observed a flare that went 
up to ~1 X lO'^^ erg s~^ and followed a roughly expo- 
nential decay from the initial detection with an e-folding 
time of ^^3 days (judging from our limited sampling). 
We note that the luminosity can very possibly have be- 
come higher, since the true peak of the flare could easily 
have been missed due to the rather short duration of the 
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event and the limited sampling. The duration of the flare 
was at least '^lO days and may have been ~10-20 days 
longer, but it cannot have started more than two weeks 
before the initial detection. The first observation of our 
monitoring program in 2010 April also showed a clearly 
elevated luminosity of ^2x lO"^"* erg s~^; this observation 
may have been made during the decay of a similar flare, 
although this is clearly highly speculative. 

The most natural explanation for these increases in 
flux is that they are accretion events. Activity at lu- 
minosities <10'^^ erg s~^ is for most Galactic transients 
hard or impossible to detect with all-sky monitors, and 
dense monitoring of quiescent sources with pointed ob- 
servations is rare. Our knowledge is therefore limited 
on how common such low-level activity is in transient 
sources and what the typical properties of such activ- 
ity are. A few tens of low-luminosity transients, often 
referred to as very faint X-ray transients (VFXTs; usu- 
ally defined as having peak 2-10 keV X-ray luminosities 
in the range ^lO^'^ -lO'^^ erg s~^), have been identified 
in the Galaxy (e.g.. iWiinands et al.ll2006t iHeinke et al.l 
I2009at IDegenaar fc Wiinand£ll2009ll2010D . A significant 
fraction of these have shown type I X-ray bursts and 
are therefore accreting neutron stars, most likely with 
low-mass companio ns (IDegenaar fc Wiinandi 120101 and 
references therein). IDegenaar fc Wiinandsl (|2009l I2010D 
report on a 4 yr monitoring campaign of the Galactic cen- 
ter with almost daily Swift observations for ^^9 months 
of the year. They detected eight faint transients show- 
ing a variety of low-level activity (some showed peak 
luminosities above the VFXT range). Four of these 
sources showed flares similar to the one we observed 
from XTE J1701-462, with durations of ~l-2 weeks and 
peak luminosities in the range ~(0.7-2) x 10^^ erg s~^; 
we note that for all four sources only one such short 
flare was observed during ~4 yr of quasi-daily monitor- 
ing, whereas our observations suggest that such events 
are more frequent in XTE J1701— 462 . We also note 
that IDegenaar fc Wiinandsl (|2009ll2010[ ) quote luminosi- 
ties in the 2-10 keV band; the peak 2-10 keV lumi- 
nosity we observed for the XTE J1701— 462 flare — with 
much sparser sampling, however — is ~5 x 10^"* erg s~^. 
It has been pointed out that the low long-term average 
mass accretion rates implied by observations of VFXTs 
may pose difficulties for binary evolution models in ex- 
plaining th eir existence (e.g.. King fc Wiinand"i 120061 : 
IDegenaar fc Wijnands 2009, 2010).' The observed be- 
havior of XTE J1701— 462 supports one simple scenario 
which could possibly explain this for some sources: if 
they undergo large outbursts with long recurrence times 
(perhaps decades or longer) and with sporadic low-level 
activity in between, then their actual long-term average 
accretion rates could be much larger than implied by a 
few years of monitoring. However, this is unlikely to 
be the case for all such sources, given that none of the 
known VFXTs have been observed to undergo a large 
outburst. We also note, as has been pointed out (e.g., 
IDegenaar fc Wiinand£ll2010l ) , that activity with such low 
peak luminosities presents a challenge to the disk insta- 
bility model th ought to describe acc r etion cycles in tran- 
sient LMXBs (|King fc Ritteii 119981 : iLasotal I2001D . and 
it is not clear that the same mechanism is at work in 
the flares from XTE J1701— 462 as in the low-luminosity 
activity of the VFXTs (nor is it clear whether a single 



mechanism is at work in all VFXTs). 

As discussed above and in Paper I, the increase in lu- 
minosity seen in our third XMM-Newton observation re- 
sulted from an increase in both thermal and nonther- 
mal flux (and this is probably also the case for the 
elevations in flux seen with Swift). The physical ori- 
gin of the nonthermal flux or the extra thermal flux is 
unclear. A prominent hard (power-law) component is 
seen in spectra from NS-LMXBs in their low-luminosity 
hard state; this component is often attributed to some 
sort of Comptonization (e.g., IDi Salvo fc Stellal [2002]; 
[Done et al. 2007). The power-law component we see dur- 
ing flaring may be a lower-luminosity version of the hard 
state component. It is also unclear whether the vari- 
able nonthermal flux seen in all the Chandra and XMM- 
Newton spectra has the same origin as the much stronger 
nonthermal flux seen in XMM-3 and the brighter Swift 
observations. The extra thermal flux observed during 
flaring may be produced as accreted matter hits the neu- 
tron star surface, perhaps in some sort of boundary layer 
between the surface and an accretion disk (similar to 
what is seen i n NS-LMXBs at higher luminosities). We 
also note that IZampieri et al.l (|T995) show low-level ac- 
cretion onto a neutron star surface can under certain 
conditions result in a spectrum with a thermal (hard- 
e ned blackbody-like) sh ape. 

iCackett et al.l (|201 01j) show that the thermal flux from 
the quiescent NS-LMXB Cen X-4 has varied somewhat 
irregularly over a period of ~15 yr; the source also has 
a variable nonthermal component. The origin of this 
variability is unclear, but the authors argue that the 
most likely explanation is that the variability in both 
the thermal and nonthermal components is due to vari- 
able low-level accretion. This conclusion is supported by 
the fact that the two components seem linked together, 
with their flux ratios staying appr oximately constant. Al l 
the quiescent spectra analyzed by ICackett et al.l ()2010bD 
show a thermal flux fraction of ~50%-60%. This is sim- 
ilar to what we observed in XMM-3, which showed a 
nonthermal fraction of ~53%, but that may well be a 
coincidence. The other quiescent observations of XTE 
J1701— 462 have shown power-law fractions ranging from 
a few percent to ~50%. As discussed in Section [521 the 
limited number of counts in our Swift spectra do not al- 
low us to place strong constraints on the separate behav- 
ior of the thermal and nonthermal components during the 
flaring we observed. Our simple hardness ratio analysis 
indicates that some of the increase can be attributed to 
the thermal component (in agreement with XMM-3), but 
that at the highest count rates observed the power-law 
component contributes the majority of the flux. We also 
note that Cackett et al. (2011), analyzing archival obser- 
vations of Aql X-1 during quiescence, found a flare with 
a likely duration of ^60 days and a ^5-fold increase in 
flux; the data suggest (but are not conclusive) that both 
the thermal and nonthermal flux rise during the flare. 

3.3. Effects of Low-level Activity on the Equilibrium 
Temperature 

Ahhough XTE J1701-462 may stiU be coohng, the 
equilibrium surface temperature is probably not very 
much lower than the current one, given how slow the 
possible cooling has become at this point. This suggests 
that the equilibrium thermal luminosity of the source is 
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high compared to other ne utron star transients (see, e.g., 
iHeinke et al.ll2007ll2009b[ ): by extension, the equihbrium 
crust and core temperature would then be comparatively 
high, unless the crust microphysics is cons iderably differ- 
ent from that used i n calculations to date (Shterni n et al.l 
[2007t [B rown fc Gumming 2009,) . It is natural to ask 
whether low-level accretion between large outbursts can 
to some extent explain the high temperature. To gauge 
the plausibility of this, we compare an estimate of the 
long-term average luminosity due to low-level activity 
to that of bright outbursts; this is equivalent to compar- 
ing the long-term average mass accretion rates, assuming 
that the radiative efficiency of the accretion is the same. 
To estimate the fluence of the Swift flare, we make the 
conservative assumptions that we caught the flare near 
its peak and that the pre- and post-peak contributions 
to the total fluence are of similar size. We also assume 
that the decay of the flare followed the best-fit exponen- 
tial mentioned in Section [2.2l Furthermore, we make the 
(admittedly highly uncertain) assumption that the typi- 
cal recurrence time for such fiares is 4 months, based on 
our speculation that the first Swift observation was made 
during the decay of a similar fiare. To m ake a rough bolo- 
metric correction, we use the result of lin 't Zand et all 
(|2007D . who in the literature on broadband spectra of 
LMXBs find a typical factor of 2.9± 1.4 between reported 
0.1-100 keV unabsorbed fluxes and 2-10 keV absorbed 
fluxes. For our three brightest Swift observations and for 
XMM-3, we flnd ratios between the unabsorbed fluxes 
in the 0.5-10 and 2-10 keV bands in the range ~2.0-2.2 
(using absorbed 2-10 keV fluxes gives a range of ~2.4- 
2.7). Since for XTE J1701-462 the absorption is high, 
we base our correction on the unabsorbed 2-10 keV flux 
values we measured and use a correction factor of 1.4 
for our 0.5-10 keV luminosities; this number is obviously 
subject to a large uncertainty, but our overall conclu- 
sions are not very sensitive to this. After subtracting 
the baseline bolometric thermal emission (based on the 
most recent Chandra observation), we find a long-term 
average fiare luminosity of ^2 x 10^^ erg yr~^. Further- 
more, assuming that the nonthermal fiux we have ob- 
served throughout quiescence is due to residual accre- 
tion (which may not be the case), and taking an av- 
erage of the fiux values derived from our Chandra and 
XMM-Newton observations (excluding XMM-3), we get 
an average luminosity of ~7 x 10*° erg yr~^, i.e., roughly 
a third of the estimated flare contribution. This then 
gives a total of ~3 x 10*^ erg yr~^ due to low-level ac- 
tivity. The total bolometric energy output of the 2006- 
2007 outburst is estimated to have been '--^1 x 10*^ erg 
(see Paper I). For the low- level activity to be a signif- 
icant factor in the total thermal budget of the neutron 
star would therefore require a recurrence time of at least 
~10* yr for large outbursts like the 2006-2007 one. A 
recurrence time this long is implausible; assuming that a 
mass accretion rate M results in a bolometric luminosity 
of Lboi = iMc^, with e = 0.2, this would imply a long- 
term average mass accretion rate in the system of a few 
times 10~^^ Mq yr~^, which is hard to reconcile with the 
high observed surface temperature (see the discussion in 
Paper I on the range of recurrence times consistent with 
standard cooling in the neutron star). 
It is interesting to note that of the four transients 



studied bv lDegenaar fc Wiinandl (|2009ll2010D that have 
shown short flares similar to the one in XTE J1701— 462, 
two have also shown longer and more intense activ- 
ity For example, GRS 1741.9-2853, a conflrmed NS- 
LMXB, showed both a ~1 week flare with estimated 
peak/average 2-10 keV luminosities of ~7 x 10'^'*/3 x 
IQ^i gpg g-i g^jjj week outburst with peak/average 

luminosities of ^2 x 10^^/1 x 10'^^ ergs"^, as well as 
an even brighter 4-5 week outburst with luminosities 
of ~1 x 103'^/2 X 10^*^ erg s'^ (we note again, however, 
that the '^l week flare is the only such very short and 
faint fiare that has been detected from the source during 
four years of monitoring). Based on the detection history 
of the source in the past decade, iDegenaar fc Wiinandsl 
(|2OT0h estimate a recurrence time (not including the 
~1 week fiare) of ~2 yr. If XTE J1701-462 experi- 
ences similar small outbursts every two years or so with 
a duration of a few months and an average luminos- 
ity of ^10"^^ erg s~^ (any past activity with luminosities 
<5 X 10"^^ erg s~^ would have been hard or impossible to 
detect with the RXTE All-Sky Monitor), then a recur- 
rence time of several hundred years for large outbursts 
would allow these smaller ones to contribute significantly 
to the average mass accretion rate. A recurrence time of 
this order is long but not inconceivable; this would im- 
ply a total long-term average mass accretion rate between 
10^^* and 10^^° M© yr~*. However, low-level accretion 
is unlikely to have produced a very large increase in the 
equilibrium surface temperature, since the fourth-power 
dependence of the bolometric thermal luminosity on the 
temperature implies that even if the average low-level 
accretion rate were of similar size to that from large out- 
bursts, this could not result in more than a ^20% higher 
equilibrium temperature than in the absence of any low- 
level activity. We therefore conclude that it is implausi- 
ble that low-level accretion can, for example, explain the 
likely '^2 times higher equilibrium surface temperature 
of XTE J1701-46 2 compared to those of MXB 1659-29 
and KS 1731-260 (jCackett et al.l[2?M[20T0al) . This dif- 
ference is more likely mostly due to a higher long-term 
average accretion rate of long-duration outbursts in XTE 
J 170 1—462 and/or more efficient cooling in the neutron 
stars in MXB 1659-29 and KS 1731-260. 

3.4. Effects of Low-level Activity on the Cooling Curve 

In our discussion, we have assumed that the thermal 
component in the nonflare quiescent spectra arises en- 
tirely from thermal emission from the neutron star sur- 
face due to heat from the crust. An important issue that 
needs to be addressed is whether the thermal component 
outside flares can to some extent be due to ongoing low- 
level accretion. A signiflcant contribution from residual 
accretion to the thermal flux observed would mean that 
the surface temperatures inferred would not accurately 
reflect the thermal state of the crust. We see increased 
thermal emission (along with a rise in nonthermal emis- 
sion) during flares which are almost certainly accretion 
events; since we see some nonthermal emission during 
regular quiescent observations as well, one might sus- 
pect that some of the thermal emission observed in those 
observations might stem from accretion. However, two 
characteristics of the observed emission indicate that this 
is unlikely to be a large effect. First, the evolution of 
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the inferred temperature throughout quiescence (exclud- 
ing the flares) is weh behaved, having shown a smooth 
monotonic decrease (with the exception of CXO-4 and 
possibly CXO-7, although the latter is consistent with 
being a statistical fluctuation), with little variation in 
thermal flux over the last seven observations (^30%). If 
accretion were in general a significant or even dominant 
contributor to the thermal emission, much more irreg- 
ular variability would seem likely; this is what we have 
observed from the nonthermal emission, which has varied 
by a factor of '^10 throughout the nonflare observations 
and by a factor of ~6 in the last seven observations. 
Second, there does not seem to be any correlation be- 
tween the thermal and nonthermal flux observed outside 
the flares, whereas both components rise together dur- 
ing fiares. It is hard to see why this behavior would be 
different outside the flares (i.e., at, presumably, lower 
accretion rates) compared to during flares if accretion 
causes a significant part of the thermal emission in both 
cases. Overall, we conclude that accretion probably does 
not play a significant role in the nonflare thermal emis- 
sion. However, we acknowledge that this is a possibility 
that cannot be conclusively ruled out and could affect 
our results to some extent. 

Another important (and separate) question is whether 
low-level activity can have a significant effect on the ther- 
mal state of the crust and thereby on the thermal surface 
emission. We first consider the effect — arising from deep 
crustal heating — of the persistent (but variable) low- level 
accretion possibly implied by the nonthermal compo- 
nent that has been observed throughout quiescence. We 
convert the observed average power-law flux to a mass 
accretion rate as described above, and assume a total 
heat deposit from nuclear reactions in the neutron star 
crust of :^ 1.9 MeV per ac creted nucleon, as calculated by 
iHaensel fc Zdunil3 (|2008| ) for a heavy-element envelope 
on the neutron star composed of ^^Fe nuclear burning 
ashes. Part of the energy produced will escape the star 
via neutrino emission, but we can get an upper limit on 
the photon emission by assuming that all of the energy is 
radiated away at the neutron star surface. This gives a 
thermal luminosity contribution of ~2 x 10^^ erg — 
compared to ^5 x 10^'^ erg s~^ currently observed — and 
corresponds to a completely negligible temperature in- 
crease of <0.1 eV. We emphasize that the thermal emis- 
sion we are considering here is due to deep crustal heat- 
ing and is not thermal emission instantaneously produced 
during accretion. Any heating of the neutron star at the 
surface due to accretion should be shallow in extent and 
thus radiated away very quickly, and would not be ex- 
pected to have long-term effects on the temperature of 
the crust. 

As estimated in Section 13. 3[ the long-term average 
mass accretion rate due to flares of the sort we observed 
with Swift should only be a factor of ^^3 higher than 
that from the possible persistent low-level accretion (cal- 
culated above), and should therefore on average have a 
very small effect on the crust temperature. However, in 
this case, the accretion comes in shorter and more inten- 
sive bursts, which could lead to temporary but larger in- 
creases in thermal flux from the surface. Calculating the 
effect of short-duration accretion events on the surface 
temperature is complicated by many factors. The nu- 



clear reactions which contribute to the total heat deposit 
per nucleon take place at various depths in the crust and 
the thermal diffusion time to the s urface varies greatly, 
from days to many hundreds of days (jBrown fc Cummingl 
l2009f l. The strongest heat sources, producing the major- 
ity of the tot al heat deposit, are at densities of 10^^- 
10" g cm-3 (jHaensel fc ZdunikI [20081 ). at depths from 
which the thermal diff usion time to the surf ace is likely 
hundreds of days (Bro wn fc Cummingl l2009h . Diffusion 
will also cause the effect of any individual reaction on the 
surface temperature to be spread out in time. In addi- 
tion, part of the heat will flow into the core rather than to 
the surface; some fraction will be radiated away through 
neutrinos. Making the same assumptions as before about 
the fluence of the flare we observed, we estimate the to- 
tal mass accreted during the flare to be ~4 x 10^° g; this 
gives a total heat deposit of ^7 x 10"^^ erg in the crust 
if we again assume that each nucleon contributes ~1.9 
MeV. We can make a very crude order-of-magnitude es- 
timate of the upper limit to the luminosity increase due 
to the flare by assuming that all this heat is radiated 
at the surface on the typica l thermal diffusion timescale 
of the crust, ~10^ days (Bro wn fc Cu mming 2009). This 
gives a luminosity of ^10'^^ erg s""'^, which would imply a 
temperature increase of ^0.5 eV — virtually undetectable 
with current instruments. We note, however, that a to- 
tal amount of accreted mass that is one or two orders 
of magnitude higher — conceivable for longer and more 
intense low-level activity — would place this upper limit 
h igh enough for a detectab le effect on the temperature. 

lUshomirskv fc Rutledgtj (|200lD study the time- 
variable quiescent luminosity of a neutron star un- 
dergoing very short (1 day) accretion outbursts with 
recurrence times of 1 or 30 yr, and find that this 
luminosity depends sensitively on the microphysics 
of the crust and core. They consider both standard 
and rapid core cooling, as well as both high and low 
thermal conductivity for the crust. In their simula- 
tions, even a total accreted outburst mass three orders 
of magnitude higher than what we estimate for the 
observed flare in XTE J1701— 462 does not result in a 
luminosity increase — due to crustal heating — of more 
than a few times 10^^ erg s~^; this holds regardless 
of what assumptions they make for the core cooling 
and crust conductivity. This suggests that periods of 
low-level activity that are considerably brighter and 
longer than what we have observed in XTE J 170 1—462 
would not be able to have an appreciable effect on the 
surface temperature. We stress, however, that there 
is considerable uncertainty in the microphysics of the 
neutron star crust and core, and we cannot exclude 
the possibility that more refined calculations — and 
calculations more tailored to the observed behavior of 
XTE J 170 1-462— would lead to different conclusions. 

The fourth Chandra observation indicates increased 
thermal flux without an associated rise in nonthermal 
flux, and the preceding observation (XMM-3) strongly 
indicates ongoing accretion. It is therefore natural to 
ask whether the possibly elevated surface temperature 
observed in CXO-4 can be attributed to crustal heating 
from the presumed XMM-3 accretion event. We conclude 
that this is unlikely given the results discussed above 
and how close in time XMM-3 and CXO-4 were. These 
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observations were only separated by ~73 days and any 
enhanced activity cannot have started more than ^^120 
days before CXO-4, whereas most of the crustal heating 
takes place at depths from which the thermal diffusion 
time to the surface is thought to be hundreds of days, as 
mentioned above. 

4. SUMMARY 

We have presented a new Chandra observation, made 
in 2010 October, of the transient NS-LMXB XTE 
J1701— 462, which entered quiescence in 2007 August af- 
ter an extraordinarily luminous 19 month outburst. This 
observation extends our monitoring of the source in qui- 
escence from ~800 days to ~1160 days since the end 
of the outburst. We have also presented the results of 
a Swift monitoring program of XTE J1701— 462, which 
took place during 2010 April-October, and whose pur- 
pose was to investigate possible low-level activity from 
the source. 

The new Chandra observation indicates that the effec- 
tive surface temperature of the neutron star may have 
decreased since the preceding Chandra observation(s) in 
2009, implying that the neutron star crust may still be 
slowly cooling toward thermal equilibrium with the neu- 
tron star core after having been heated during the 2006- 
2007 outburst. An additional observation further into 
quiescence is needed to conclusively determine whether 
cooling is still ongo i ng, an d if so, to constrain its rate. 
iBrown fc Cummind ()2009[ ) calculate that the form of 
the cooling curve should be a broken power law level- 
ing off to a constant at late times. With the present 
data, the shape of the overall cooling curve is consistent 
with a broken power law. However, the observed break 
(at ~25-80 days into quiescence) is much earlier than 
the break predicted by theory (at a few hundred days 
post-outburst) and may therefore have a different origin. 
We can also not exclude that the cooling has followed 
an exponential decay to a constant level, although on- 
going cooling would make this somewhat unlikely. The 
rapid initial cooling during the first ~200 days of qui- 
escence strongly indicates a highly conductive neutron 
star crust, and possibly suggests low-impurity m aterial 
(jShternin et al.l 120071 IBrown fc Cummind [20091 ). The 
high current surface temperature (which corresponds to 
a bolometric thermal luminosity of ^5 x 10'^^ erg s~^), 
coupled with the slow (or possibly nonexistent) current 
cooling, suggests that the equilibrium crust / core temper- 
ature of XTE J1701— 462 is high compared to other tran- 
sient NS-LMXBs, unless the microphysics of the crust is 
consi derably different from that used in calculation s to 
date (jShternin et al.ll2007l: [Frown fc Cummindl2009D . 

We observed a large temporary increase in luminos- 
ity with XMM-Newton ~230 days into quiescence. This 
prompted us to undertake a 7 month Swift monitoring 
program of XTE J1701— 462 with short observations once 
every two weeks. During this time, we observed a short- 
term (^^10-20 day) flare — presumably arising from low- 
level accretion — which went up to a luminosity of at least 
~1 X 10^^ erg s"""^, i.e., ^20 times higher than the nor- 
mal quiescent level. We also observed a smaller increase 
in luminosity at the beginning of the program. Accre- 
tion during events like these will lead to some heating in 
the crust. A simple analysis suggests that flares of the 
magnitude observed are not likely to have significantly 



affected the equilibrium crust/core temperature of the 
neutron star and should not be able to have a signifi- 
cant impact on the cooling curve. However, it is possible 
that brighter and longer periods of low-level activity — as 
have been observed in some faint Galactic neutron star 
transients — may have had an appreciable effect on the 
equilibrium temperature of the neutron star. It is im- 
portant that the source be monitored further to better 
constrain its low-level activity in quiescence. 
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